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ABSTRACT 

We obtained 3 epochs of simuhaneous Hubble Space Telescope (HST) /Wide Field Camera 3 and 
Chandra observations of the ultraluminous X-ray source NGC 5408 X-1. The counterpart of the 
X-ray source is seen in all HST filters, from the UV through the NIR, and for the first time, we 
resolve the optical nebula around the ULX. We identified a small OB association near the ULX that 
may be the birthplace of the system. The stellar association is young, ~ 5 Myr, contains massive 
stars up to 40 Mq, and is thus similar to associations seen near other ULXs, albeit younger. The 
UV/Optical/NIR spectral energy distribution (SED) of the ULX counterpart is consistent with that 
of a BOI supergiant star. We are also able to fit the whole SED from the X-rays to the NIR with an 
irradiated disk model. The three epochs of data show only marginal variability and thus, we cannot 
firmly conclude on the nature of the optical emission. 

Subject headings: accretion, accretion disks - black hole physics - X-rays: binaries - X-rays: individual 
(NGC 5408 X-1) 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are extragalactic 
sources that are not at the nucleus of their galaxy and 
emit well above the Eddington limit of a 20 Mq black 
hole {Lx 3 X lO'^^ erg s""'^). The most important, still 
unresolved question is whether these objects contain in- 
termediate mass black holes - IMBHs - (M > 100 Mq) or 
stellar mass b lack holes with supe r-Eddington or beamed 
emission (see lFeng fc Soriall2011l for a review). 

A direct answer regarding the mass of compact objects 
in ULXs would be possible if a radial velocity curve of 
absorption lines from the companion star, or emission 
lines coming from the accretion disk, could be obtained. 
To date, this has been proven to be a real challenge due 
to the lack of absorption lines in ULX optical spectra and 
what seem to be non-periodic variations of the only trace- 
able emission lin e. He IIA4686, in Ho lmberg IX X-1 and 
NGC 1313 X-2 (iRoberts et al.|[20lTh . The only excep- 
tion, to date, is the optically ver y luminous source P13 
in NGC 7793 (jMotch et al.ll2011| ) where absorption hues 
from a late B supergiant star (in addition to a He II emis- 
sion line) have been observed. This is a very promising 
candidate to constrain the mass function of the system. 
However, this is the only example out of more than ten 
identified optical counterparts. 

In all other ULXs, the true nature of the donor star is 
mostly unknown. The X-ray lightcurve of M82 X-1 shows 
a 62 day periodicity consistent with a strictly periodic 
signal a nd hence is likely the orbital period of the binary 
system (jKaaret fc Fenej|2007[ ) . If the companion star fills 
its Roche lobe, this means that such star has to be of low- 
density and hence would be a giant or supergiant star 
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(jKaaret. Simet. fc Lang|[2006[ ). Other constraints on the 
nature of the donor stars rely on studies of their envi- 
ronment. Indeed, some counterparts belong to or arc lo- 
cated nearby small steU ar associations (Soria et al. 20051; 
IGrise et al.lf2008L [20TTh . where it is possible to constrain 
the association age, and thus the maximum mass for a 
member star. This has led to maximum mass estimates 
of 15 — 20 M0 that apply as well to the donor star of 
the ULX systems, if mass transfe r did not play an im- 
portant role. In some other cases (|Feng fc Kaaretl[200a 
I Yang et al.ll2011| ). the constraints arise from field stars in 
the absence of such associations. But these estimates are 
based on indirect evidence and the precise stellar type of 
those stars remains unknown. 

The most obvious problem in identifying the compan- 
ion star in the optical emission of ULXs is that the 
flux contribution coming from the accretion disk (direct 
and irradiated) is not known but is definitely present in 
some sources. The best hint comes from NGC 1313 X-2, 
where its optical lightcurve is contaminated by short- 
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Broad He II emission lines with equivalent widths up 
to ten times higher than in the brightest Galactic high- 
mass X-ray binaries are present in the optical spect rum of 
some counterparts (e.g . iPakull et al. 2006: Robert s et al.l 
[201ltlCseh etaII[20Tiri and are also a direct proof of sig- 
nificant X-ray reprocessing in these systems. 

The counterpart optical colors have often been used 
to place constraints on the companion star. The high 
optical/UV luminosi ties and blue colors suggest early - 
type, OB stars ( e.g. iLiu et al.ll2002l : iKaaret et al."2004 
iSoria et al.l 120051) . But as suggested bv IKaaret. (200i 
and IPakull et all ([200l . the optical light may be con- 
taminate d, or even dom inated, b y the irradiated accre- 
tion disk (jKaaret fc Corb el 20091: " lTao et~aI]|201lD . This 
is similar to active low-mass X-ray binaries, although 
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ULXs have higher X-ray luminosities and may have more 
massive donor stars. Unfortunately, the optical spectral 
energy distributions of hot, irradiated disks look rather 
like those of OB stars. 

Our present work is devoted to a well studied ULX, 
NGC 5408 X-1 (hereafter N5 408X1), located at 4.8 
Mpc (jKarachentsev et all |200^ . N5408X1 is a bona- 
fide ULX with a n average luminos ity (0.3-10 keV) of 
1 X 10'*° ergs-i (|StrohmaveJ 120091) that displavs clear 
variability within a factor of ^ 2 |Kaaret &: Fengj 120091 : 
lStrohmaveill2009t lKongj|20lH) . Intensive S'm/f monitor- 
ing of the source has been used to infer a ^ 115 day pe- 
riod that has been interpreted as the orbital period of the 
ULX system (|Strohmavedl20090 . although [Foster et all 
(|2010( ) suggests that this modulation may instead be su- 
perorbital. The X- ray spectrum of N540 8X1 docs not 
seem to vary much (jKaaret fc Fend [2009( 1 . All observa- 
tions performed to date show a soft spectrum F 2.5 and 
the presence of a soft component. Different interpreta- 
tions have been put forward: if interpreted in the frame- 
work of Galactic black holes binaries (GBHBs), stan- 
dard models (powerlaw plus disk component) point to 
the presence of a massive black hole if the t emperature of 
the d i sk scales as i n Galactic counterparts (iKaaret et al.l 
[200l ISoria et al.l l200l IStrohmaver et al.l I2007D . The 
power-law would be in that case the signature of a corona 
as seen in GBHBs. Instead, the interpretation in the 
context of super-critical accreting systems would mean 
that the black hole is much smaller, with the possibil- 
ity that a strong wind is the origin of the soft compo- 
nent, while the high energy part of the spectrum would 
be due to disk emissi on inside the spherization radius 
(jPoutanen et al.ll2007D . Alternatively, the soft compo- 
nent may represent emission from the outer disk with the 
hard emission associated with thermal Comp tonization 
(jGladstone et al.|[200l [Middleton et al.ll201lD . Finally, 
the soft excess may arise from re flection from the disk 
(|Caballero-Garcfa fc Fabia^lMol l. 

N5408X1 is also one of the few ULXs where quasi- 
periodic oscillations (QPOs) have bee n discovered 
i|Strohmaver et al.ll2007l: lStrohmave^l2009D . By analogy 
with Galactic black hole binaries (GBHBs). IStrohmaved 
(|2009f l have interpreted this feature as a type C low- 
frequency QPO (LFQPO), w hich would i n dicate that 
N5408X1 hosts an IMBH. Mi ddleton et al.l (|2011| ) have 
challenged this conclusion and concluded that even the 
most super-Eddington BHB, GRS 1915+105, is not re- 
ally a good match to the properties of N5408X1. They 
hypothesized that the temporal and spectral properties 
are more consistent with the feature seen in some Nar- 
row Line Seyfert 1 galaxies which are supposed to be even 
higher super-Eddington sources, which would indicate a 
black hole with mass < 100 Mq. But the large uncer- 
tainties associated with the black hole masses in those 
systems avoids any definitive conclusion. 

N5408X-1 is also seen at other wavelengths. This 
source is one of the few ULXs detected in radio, the 
emission being con sistent with a powerful nebula sur- 
rounding th e UL X (iKaaret et al.ll200l ISoria et al.ll2006l: 
iLang et al.ll2007t iCseh et al.ll201lD~This nebula is also 
detected in optical wavelengths d Pakull fc Mirionil 120031 : 
ISoria et al.l 120061 IKaaret fc CoTb el 2009) where the op- 
tical emission is consistent with X-ray pliotoionization. 



The optical counterpart of the ULX, first uniquely iden- 
tifi ed in a Hub b le Spa ce Telescope (HST) /WFPC2 image 
bv iLang et al.l ()2007D as a Vb = 22.2 magnitude object, 
ha s been subsequent l y stu died using VLT s pectro scopy 
bv IKaaret fc Corbell ((20091) and ICseh et all (|20TTI) . The 
optical spectrum of N5408X1 is similar to many other 
ULXs, displaying a blue continuum with superimposed 
strong and narrow emission lines coming from the neb- 
ula. Of prime interest is the presence of broad compo- 
nents with widths of ~ 750 km s~^ in the He II and Hp 
lines (jCseh et al.l 1201 ll ) that appear spatially point-like, 
similar to the optical continuum of the counterpart. The 
radial velocity of the broad He II line varies, which may 
represent the motion of the compact object or m otions 
within the accretion disk (see iRoberts et al.|[20Tll for the 
behavior of this line in other sources). 

In this paper, we present new, simultaneous Chan- 
dra/ HST observations of N5408X1. HST observations 
that allow us to present for the first time the spectral 
energy distribution of a ULX from the ultraviolet (UV) 
to the near-infrared (NIR) wavelengths. Combined with 
simultaneous X-ray data, we are able to show that an 
irradiated disk model fits well all the data. We discuss 
other alternatives and suggest a way to confirm this re- 
sult with new observations. 

2. OBSERVATIONS AND DATA ANALYSIS 

We obtained three series of multiwavelength observa- 
tions using the recently installed Wide Field Camera 3 
(WFC3) instrument aboard the Hubble Space Telescope 
(HST) (Program ID 12010, PI: P. Kaaret, along with 
simultaneous Chandra X-ray observations (Program ID 
11400085, PI: P. Kaaret). 

2.1. Optical observations 

Each HST visit consisted of observations in UV, op- 
tical, and near-infrared filters. UV and optical observa- 
tions were carried out in the F225W, F336W, F547VI, 
and F845M filters using the WFC3/UVIS camera. We 
used the WFC3/IR camera for the near-infrared obser- 
vations, in the F105W and F160W filters. A single ob- 
servation was done with the WFC3/UVIS using the nar- 
row F502N filter to isolate emission due to the forbidden 
[O III]A5007 line. A summary of the observations can be 
found in Table [J 

We performed aperture and point-spread function 
(PSF) photometry on the drizzled images. Aperture 
photometry was done only on the four UVIS images, 
since the crowding is really severe in the infrared images, 
preventing any reliable photometry in the region nearby 
the ULX with this method. We perfor med the aperture 
photo metry with Sextractor v. 2.5.0 CBcrtin fc Arnout^ 
11996( 1. using an aperture of 2 pixels in radius (i.e, 0.08"). 
The background was estimated by relying on the global 
mapping done by Sextractor. Aperture corrections were 
calculated using relatively bright and isolated stars, per- 
forming photometry using a 0.4" radius (i.e. 10 pixels). 

PSF photom etry was performed using DAOPHOT II 
(|Stetsonlll987[ l under ESO-MIDAS. For the UVIS images, 
we used a 2.5 pixel aperture (i.e. 0.10") which corre- 
sponds approximatel y to 60-70% of the enclosed energy 
((Kalirai et al.ll2009bl ) depending on the filter. The back- 
ground was chosen to be a 10 pixel width annulus located 
at 15 pixels from the object centroid. Given the use of 



The ultraluminous X-ray source NGC 5408 X-1 



3 



a sub-array for the UVIS observations, we used a con- 
stant PSF across the frame. Aperture corrections were 
calculated using the aperture photometry routine within 
the DAOPHOT package in the same way as previously de- 
scribed. A similar procedure was applied to the IR im- 
ages, with a 2.0 pixel aperture (i.e. 0.26"), correspondin g 
to 70-80% of the enclosed energy ()Kalirai et al.ll2009aD . 
The IR observations were performed in full frame mode, 
thus we used a quadratically varying PSF in the field. 
Aperture corrections were measured as for the UVIS im- 
ages. Finall y, zeropoints calculate d for a 0.4" radius were 
taken out of iKalirai et"all ()2009bl ial) for UVIS and IR ob- 
servations. 

It turns out that the main uncertainty in the photom- 
etry is due to the aperture correction. Indeed, the use 
of a window for the UVIS observations means that the 
number of bright stars is low which forces us to rely on 
fewer stars, or less bright stars, for the correction. 

In the rest of this paper, we will use the results coming 
from the PSF photometry, but the two methods give con- 
sistent results for bright or isolated stars, with differences 
below 0.1 mag. 

2.2. X-ray observations 

N5408X1 was observed with Chandra using the Ad- 
vanced CCD Imaging Spectrometer (ACIS) instrument. 
For each series of T/^T observations, we obtained a simul- 
taneous 12 ks Chandra exposure. Given the high count 
rate expected of N5408X1, only one chip (ACIS S3) was 
used, with a 1/8 subarray mode. This reduces the frame 
time to 0.441 s which mitigates pile-up in the data. The 
observed count rate in the three observations varies from 
~ 0.3 to ~ 0.4 counts/s. Standard extr action of spectra 
was d one using the CIAO 4.3 package (jFruscione et al.l 
120061 ) with CALDB 4.4.1. We rebinned the spectra to a 
minimum of a 5ct significance per bin. All fitting was p er- 
formed using ISIS v. 1.6. 1-26 (|Houck fc Denicolall2000l) . 

3. RESULTS 

The opti c al cou nterpart of the ULX, as identified by 
ILang et all (j2007D . is visible in aU the HST Sltcis used in 
this study. Figure [T] shows a composite image including 
the F225W, F502N ([O III]) and F845M fihers. Magni- 
tudes and corresponding fluxes of the ULX counterpart 
are presented in Table [2J The image of the field around 
N5408X1 (Figure [1]) highlights the presence of a small 
stellar association, close to the ULX at about 4"(~ 93 pc) 
north-east with ^ 20 bright, resolved stars. In the field of 
view of our WFC3 observations, the other nearest asso- 
ciation is at ~ 8" north-west from the ULX, and is likely 
connected to the super star clusters that are located 
norther, as s een on previous broadband HST and SUB- 
ARU images (jKaaret et al.|[2003HSoria et al.|[200l . This 
image also shows the nebula to be well resolved, with a 
diameter ~ 2. 5" (~ 60 pc) , consistent with the siz e of the 
radio nebula (|Lang et all 120071: iCseh et al.l 1201 1[) . From 
optical long-slit spectroscopy, iKaaret fc Corbell (j2OO90 
derived a size of 30 pc (full width at half-maximum) for 
the He III region along the slit. The size of the nebula in 
lower ionization emission lines, e.g., [O III], is broader, 
about 40 pc (full width at half-maximum) with broad 
wings. This is quite consistent with the HST imaging. 

3.1. Stellar environment 




Figure 1. HST/WFCS composite image (blue - F225W, green - 
F502N, red - F845M) of the field surrounding NGC 5408 X-1. 1" 
represents 23.3 pc at the distance of NGC 5408 X-1. Overplotted 
are the radio positions from VL A (yellow full cir cle) and ATCA 
observations (red dashed circle) I ILang et al.ll2007l ) with a conser- 
vative error circle of 0.4" in radius. The white box (5.5" X 4.0") 
shows the nearby stellar association. Stars located inside this rect- 
angle were plotted differently in the color-magnitude diagrams (see 
Figures [2] &|3ll. 



We used color-magnitude diagrams to investigate the 
stellar environment of the ULX. Data have been cor- 
rected by the reddening derived fro m the nebula sur- 
round ing the ULX, E{B-~V) = 08 (IKaaret fc Corbell 
I2009D using the extinction law of iCardelh et al.l (11989^ 
with i?v = 3.1. The metahicity of NGC 5408 has been 
estimated to be about a tent h of the solar metallicity, 
with 12 + log{0/H) = 7.99 (|Mendes de Oliveira etahl 
120060 . which translates i nto Z = 0.002 (base d on a stan- 
dard solar composition, iGrevesse &: Sauvall |1998() . We 
note that the age of the young and blue stars that are 
most common in the color-magnitude diagrams are not 
really sensitive to metallicity and therefore our age es- 
timate does not depend on this parameter. However, 
the stars on the red supergiant branch are more sensi- 
tive to this effect; the red supergiants apparent in the 
(F547M,F845M) diagram indicate a subsolar metallicity, 
but with a best fit consistent with Z — 0.008. 

Thus, we overplotted isochr on es from the P adova 
librarv (iBertel li et all 119941 : iMarigo et aD 120081 : 
iGirardi et al.l 12008. 201(1) with a metalhcity Z = 0.008. 
We present (Figure [2]) three diagrams using the F225W, 
F336W, F547M and F845M fihers. The diagrams with 
blue colors (upper panels of Figure[2]) show that the stars 
associated with the relatively dense association nearby 
the ULX form a young upper main-sequence/supergiant 
branch that is well fitted by an isochrone of age 
~ 5 Myr. A young age for the association is supported 
by the presence of six bright {Mf547m ~ Afy = —6 
to —7) stars with blue colors that are most likely blue 
supergiants. Again, we emphasize that the age of these 
young stars does not depend much on the metallicity, 
using Z = 0.008 instead oi Z = 0.002 introduces at 
most an uncertainty of 1 Myr which is smaller than the 
uncertainties due to photometric errors and isochrones 
models. We also overplotted Padova evolutionary tracks 
with various masses on the color-magnitude diagrams 
(Figure [3]). From those, it is apparent that numerous 
stars from the association (and from the field) are 
massive objects with masses up to 30-40 Mq. 



4 Grise et al. 



Table 1 

The HS T/WFC3 observations for NGC 5408 X-1. 





TT^ 


Instrument 


Filter 


Uatc 


Exposure time (s) 


Epoch 1 


IDQCUiUvSU 


U Vlb 


r ZZO W 


zUlU May Uz 


560 




IDQCUiUzU 


U Vlb 


r oooW 


oni n i\/Tr>,r no 
zUlU May Uz 


280 




i U QUU -L U -L U 


U V lO 


■pc:/i 7i\/T 


oni n A/To^r 09 
zuiu iviay uz 


ZUU 




ibde01040 


UVIS 


F845M 


2010 May 02 


280 




ibde01060 


IR 


F1 05 W 


2010 May 02 


298.5 




ibdeOlOSO 


IR 


FlfiOW 


2010 May 02 


498.5 


rLpocn z 


ibacUzUoU 


U Vio 


r ZZO W 


oni n TV/Tri^T ic; 
zUlU May io 


560 




iDaeUzUzU 


U Vlb 


r oodyV 


oni n TV/Tr-^T ic; 
zUlU May io 


280 






U V lo 


r O'-t i ivi 


oni n A/To^r 1 ^ 

ZUlU iVl dy _L o 


900 
ZUU 




ibde02040 


UVIS 


F845M 


2010 May 15 


280 




ibde02060 


IR 


F105W 


2010 May 15 


298.5 




ibde02050 


IR 


F160W 


2010 May 15 


498.5 


Epoch 3 


ibde53030 


UVIS 


F225W 


2010 Sep. 12 


560 




ibde53020 


UVIS 


F336W 


2010 Sep. 12 


280 




ibde53010 


UVIS 


F547M 


2010 Sep. 12 


200 




ibde53040 


UVIS 


F845M 


2010 Sep. 12 


280 




ibde53060 


IR 


F105W 


2010 Sep. 12 


298.5 




ibdc53050 


IR 


F160W 


2010 Sep. 12 


498.5 




ibde04010 


UVIS 


F502N 


2010 Dec. 26 


2600 



Table 2 

Brightness of the ULX counterpart in different filters, from the HST/WFC3 observations. Magnitudes are expressed in the HST/WFCS 
Vegamag system. The errors on the magnitude include the aperture correction uncertainty. The value between brackets is the 

photometric error, as derived by DAOPHOT. 





Filter 


Aperture correction 


VEGAmag 


Flux (10^^* erg s"^ cm"^ A"^) 


Epoch 1 


F225W 


0.44 ± 0.03 


20.09 ± 0.04 


(0.024) 


38.60 ± 1.42 




F336W 


0.42 ± 0.06 


20.64 ± 0.07 


(0.029) 


18.02 ± 1.16 




F547M 


0.53 ± 0.06 


22.43 ± 0.07 


(0.029) 


3.91 ±0.25 




F845M 


0.59 ± 0.04 


22.42 ± 0.08 


(0.067) 


1.03 ± 0.08 




F105W 


0.06 ± 0.03 


22.84 ± 0.09 


(0.072) 


0.34 ± 0.03 




F160W 


0.20 ± 0.05 


22.48 ± 0.10 


(0.091) 


0.12 ±0.01 


Epoch 2 


F225W 


0.42 ± 0.04 


20.04 ± 0.05 


(0.034) 


40.42 ± 1.86 




F336W 


0.45 ± 0.05 


20.55 ± 0.05 


(0.022) 


19.58 ±0.90 




F547M 


0.50 ± 0.06 


22.30 ± 0.07 


(0.026) 


4.41 ±0.28 




F845M 


0.54 ± 0.04 


22.37 ± 0.07 


(0.052) 


1.08 ±0.07 




F105W 


0.22 ± 0.07 


22.63 ± 0.10 


(0.090) 


0.41 ± 0.04 




F160W 


0.22 ± 0.05 


22.37 ± 0.10 


(0.073) 


0.14 ±0.01 


Epoch 3 


F225W 


0.49 ± 0.04 


19.97 ± 0.04 


(0.028) 


43.11 ± 1.59 




F336W 


0.53 ± 0.08 


20.56 ± 0.08 


(0.021) 


19.40 ± 1.43 




F547M 


0.51 ± 0.08 


22.26 ± 0.08 


(0.014) 


4.58 ± 0.34 




F845M 


0.61 ± 0.04 


22.26 ± 0.06 


(0.038) 


1.20 ±0.07 




F105W 


0.16 ± 0.07 


22.55 ± 0.12 


(0.076) 


0.44 ± 0.05 




F160W 


0.14 ± 0.06 


22.44 ± 0.10 


(0.100) 


0.13 ±0.01 



3.2. ULX 

We present (Figure 2]) a detailed spectral energy dis- 
tribution (SED) of this ULX, from 2400 to 15400 A. The 
three SEDs taken at different epochs are well consistent 
with each other. In each case, we fitted the SED by 
using a simple power- law with oc i/", F^, being the 
flux in units of erg s~^ cm^^ Hz^^. The indices of the 
power-law fits are respectively 1.37 ± 0.03, 1.33 ± 0.03, 
and 1.31 ± 0.03 for SED 1, 2 and 3. However, these fits 
are unacceptable with xl of 7.7, 7.7 and 5.0 for 4 de- 
grees of freedom. Using only the F336W, F547M and 
F845M values to constrain the fit, wc obtain a good fit 
with < 1- The indices arc not too different, with 
1.38 ± 0.08, 1.40 ± 0.07, and 1.28 ± 0.08 for SED 1, 2 



and 3 and the normalizations at A = 5500A are respec- 
tively 4.87 ± 0.16, 5.30 ± 0.15 and 5.54 ± 0.17 in units of 
10~^^ crgs~^ cm~^ . This result is c onsistent with the 
index value found bv lTao et all (1201 ID (1.41 ±0.14) from 
4 simultaneous observations obtained by HS T/WFPC2 
in April 2009 and covering the same wavelength range 
(3300-8000 A). As can be seen on Figure H those fits 
show a clear deviation from a straight power law in the 
NIR wavelengths, and also show marginal evidence for 
a turn-off at the shortest wavelengths. Interestingly, the 
deviations at these wavelengths seem to correspond to 
what is expected from a BOI star (Figure SI right panel). 
Only the flux at 15400 A does not follow the decreasing 
trend that should continue in the infrared. Wc note that 
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the power-law index of the continuum in the 4000-6000 A 
interval found by iKaaret fc COTbell (|2009l) [2.0^0 2)^ is 
likely contaminated by the nebular emission. The spa- 
tial resolution of HST/WFC3 allows us to measure the 
flux of the continuum of the counterpart, with a very low 
contamination of the nebula as can be seen in the com- 
posite image (Figure [ij that contains the nebular [O III] 
emission. 

The three X-ray spectra observed simultaneously with 
the HST data have count rates of 0.30, 0.32 and 0.42 
counts/s. Their moderate exposure times imply that 
the number of counts is limited (~ 4000) and thus 
comp lex models such as used recently in the litera- 
ture (Ipladstone et al.l 120091 : iCaballero-Garcfa fc Fabiaiil 
120101 : IWalton et al.l 120 lit ) cannot be constrained accu- 
rately. Specifically, the lack of counts at high energy 
does not allow one to se e the high-energy curvature as 
shown bv lStobbart et al.l (|2006( ). We modeled the X-ray 
spectra using basic models such as an absorbed power- 
law and a power- l aw plus disk blackbody (DISKBB, 
iMitsuda et al.l[r984f ) combination to compare with pre- 
vious results in the literature. Extinction was mod- 
eled usin g the Tuebingen-Bo ulder ISM absorption model 
fTBABS. IWilms et al.ll2000 ). We used a fixed compo- 
nent related to the Galactic extinction with nn = 0.057 x 
10^^ cm~^ (jPickev fc Lockmanlll990[) . and let a second 
extinction component varying freely. N5408X1 does not 
appear to show any strange behavior. The fits parame- 
ters for these simple models (Table Bl) are in agreement 
with previous observations (jKaiava fc PoutanenI I2009D 
where the authors argue that N5408X1 displays a L^-F 
correlation, which holds here. 

To further understand the whole SED and the ef- 
fect of X-ray emission at lower energies, we extrapo- 
lated the emission of the disk component. As shown by 
IKaaret et al.l (|2004[ ) in the case of Holmberg II X-1, some 
precautions need to be taken since the X-ray emission of 
most ULXs is largely dominated by the power-law com- 
ponent and not the accretion disk, and extrapolating the 
power law to lower energies would be non-physical. In 
the case of N5408X1, the ratio between the two compo- 
nents is about 50% ( Table IH), requiring care. Following 
iBerghea et al.l (120101 ). we modeled separately the pow- 
erlaw and the accretion disk component (Table |4]) . We 
fitted a powerlaw model to the spectra between 1 and 
7 keV, and a disk blackbody model between 0.3 and 1 
keV. As can be seen in Figure [5l the extrapolated disk 
model underestimates the HST fiuxes by about an order 
of magnitude, by a factor ~ 38/36/24 at 2400 A and 
by a factor ~ 5/5/3 at 15400 A . We also used a model 
that consists of a disk plus a thermal Comptonization 
component (DISKPN -I- CO MPTT, Table 111), as u sed in 
a recent study of N5408X1 ()Middleton et al.(l20Tll ). The 
comptonization component is a better choice since it is 
negligible in the optical wavelengths, so that there is no 
fiux contamination in the disk component, as it is the 
case for the model based on a powerlaw. Using that 
model, we found that the extrapolated disk spectrum is 
11/16/6 times lower at 2400 A and 5/8/3 times lower at 
5400 A compared to HST fiuxes. At 15400 A the extrap- 
olated disk is 1.5/2.3 lower in the two first observations 
and slightly above for the third observation compared to 
iJST fluxes (FigurejS]). We note that these extrapolations 



assume a very large outer disk radius (^ 1 x 10^"^ cm). 
Assuming narrower disks (or truncated disks) would re- 
duce their flux in the IR, optical, and UV. 

Reprocessed X-ray irradiation is likely to be significant 
in ULXs, th e optic al light curves of, e.g., N GC 1313 X-2 
(jGrise et al., .2008: ilmpiombato et all 120111) show short- 
term, stochastic variability that does not seem compati- 
ble with just the ellipsoida l variations due to acompan- 
ion star. Thus, following IKaaret fc Gorbell (|2009f ). we 
fitted the X- ray and the HST data with the irradiated 
disk model of lGierlihski et al.l ()2009( ). DISKIR. It is based 
on the standard DISKBB model but includes effects due 
to disk irradiation and Comptonization. The irradiated 
inner disk and the comptonized tail illuminate the outer 
disk, which implies a higher luminosity for this part of 
the disk. This model has b een applied successful ly on a 
galactic black hole binary (jGierliiiski et al.ll2009( ) where 
the authors conclude that the optical/UV emission in the 
soft state is consistent with reprocessing of a constant 
fraction of the bolometric X-ray luminosity. It has also 
been applied to the microquasa r GRS 1915+105 where 
the authors (jRahoui et al.l[20Tol ) show that an excess of 
mid-IR emission ('^ 4 — 8/im) is probably related to some 
reprocessed soft X-ray emission in the outer part of the 
disk. 

The DISKIR model contains nine parameters with two 
of them being constrained by the UV/optical/NIR mea- 
surements. Three parameters were frozen: the frac- 
tion of luminosity in the Compton tail which is ther- 
malized in the inner disk, was set to 0.1. The electron 
temperature and the irradiated radius were found to be 
poorly constrained so we fixed them respectively to 50 
keV and 1.1 x i?in. The electron temperature does not 
seem to have a strong impact on the other parameters. 
We tested with two other values (1.5 and 500 keV) but 
the effect on all parameters is within 20%. Changing 
the irradiated radius from 1.1 to 2.0 x i?in increases the 
inner temperature of the disk by 10%, the ratio Lc/Ld 
by 80% and fout by 20% but this does not affect our 
conclusions. The results of the fits are shown in Table 
|3| and Figure [5] and arc consistent for the three spec- 
tra. We note that spectrum 1 has worse than the 
other spectra. The residuals have positive peaks between 
0.5 an d 1 keV and could be due t o low e nergy emission 
lines. ICaballero-Garcfa fc Fabian! ()2010[ ) has suggested 
that ULXs may produce highly ionized oxygen and iron 
emission lines via refiection off the disk. However, it is 
not clear why such lines would be present in the first 
spectrum only (that has the lowest number of counts). 
This interesting issue should be looked into with deeper 
spectra. In this paper we limit ourselves to the contin- 
uum modeling. 

4. DISCUSSION 

4.1. Environment of the ULX 

ULXs are usually located inside or ne ar stellar-forming 
regions tha t look li ke OB ass ociations (jSoria et al.ll2005t 
iGrise et al..,2008 : Swartz et al. 2009 ). This has strenght- 
ened their identification with X-ray binaries contain- 
ing a massive companion star that is needed to explain 
their high, persistent X-ray luminosity. The fact that 
N5408X1 is located nearby (but not inside) to what ap- 
pears to be an OB association (Figure |T]) is thus not un- 
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Figure 2. i/ST/WFC3 color-magnitude diagrams for the stellar field (512 X 512 pixels, i.e. 20" X 20") around the ULX. Padova isochrones 
for stars of different ages are ovorplotted. Typical photometric errors are also plotted. Data have been corrected for an extinction of 
E(B — V) = 0.08 mag, the bar at the bottom right corner illustrating this effect. Upper left panel: Color-magnitude diagram in the 
(F225W,F336W) system. Upper right panel: Color-magnitude diagram in the (F336W,F547M) system. Lower left panel: Color-magnitude 
diagram in the (F547M,F845M) system. The same isochrones are plotted in the three panels, i.c for 5, 10, 20, and 50 Myr at Z = 0.008. 
Filled green circles are stars in the field, red diamonds are stars located in a rectangle of 5.5 X 4.0" coincident with a nearby OB association. 
The red star is the position of the ULX counterpart. In the bottom panel, we show as orange X's stars from the OB association that are 
in common with the upper panels. The smaller number of stars in common with the upper panels is simply due to the fact that we mainly 
see redder stars in the F845M filter that are not strong UV emitters (and vice-versa in the top panels where the strong UV emitters are 
not seen in the redder filters). 

but in any case the formation of such massive remnants 
would need rnassiy e and dense clusters. As noted by 
iKaaret et al.l (|2003[ ). the projected displacement of ~ 280 
pc from the super star clusters could be consistent with 
an ejected X-ray binary. Using r easonable ejection speed 
(10 km/s), IKaaret et all ()2003[ ) concluded that the bi- 
nary could traverse this distance in 30 million years and 
thus have moved to the present location within the life- 
time of a massive companion star. We remark here that 
30 Myr would in fact exclude most early types of stars, 
since 30 Myr is the lifetime of a ~ 10 Mq star. 

Another possibility, as noted above, is that the binary 
comes from the closer OB association located at 100 pc 



usual. For example, NGC 1313 X-2 (IGrise et al.ll2008D is 
displaced from its apparent parent cluster by 100 pc. 
This is a good argument that such X-ray binaries may 
have been ejected from their host cluster (jKaaret et alj 

W). , ^ ^ 

IKaaret et al.l ()2003l l has shown that the ULX is 
located at about 12" from the star formation re- 
gions of NGC 5408 that contain super star clus- 
ters. The runaway-co alescence scenario studied by 
iVanbeveren et al.l (l2009f ) (and references therein) is able 
to explain the formation of massive black holes up to 
several 100 M© in subsolar metallicity environments. It 
is not clear if IMBHs could be produced by this way. 
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Figure 3. HST/WFC3 color-magnitude diagrams witli Padova evolutionary tracks for stars of different initial masses (7, 10, 15, 20, 30, 
and 40 Mq) with Z = 0.008. The legend is the same as Figure [2] Upper left panel: color-magnitude diagram in the {F225W, F336W) 
system. Upper right panel: color-magnitude diagram in the (i^336iy, F5A7M) system. Lower right panel: color-magnitude diagram in the 
(F547M, F845M) system. 



from its position. In that case, it is very unhkely that an 
IMBH has been formed in such a steUar environment. We 
also note that the age of the association, ^ 5 Myr would 
imply an escape velocity of 25 — 50 km/s with the 
assumption that the binary was ejected between 1 and 
3 Myr after the stellar association was formed. These 
velocities, albeit large, cannot be ruled out since some 
OB-supergiant binaries hav e been seen with runaway ve - 
locities of that order (e.g., Ivan den Heuvel eFaLll2000l ). 
In that case, a more massive donor star would be allowed 
in the system, in comparison with the super star cluster 
hypothesis. 

4.2. UV/optical/NIR SED 

The fact that the SED of N5408 X-1 may be described 
approximately by a power-law, and does not show any 
clear clue of a peak of a blackbody leads to several con- 
clusions. As can be seen in Figure|4j only early type stars 



display a power-law extending far in the UV because the 
peak of their spectrum is located at shorter wavelengths 
than our observing limit 2250 A equivalent to a peak 
temperature of ^ 13000 K). Basically, this means that if 
the light of N5408 X-1 was dominated by a star, we could 
rule out any star of spectral type later than B5 (this is 
quite independent of the luminosity class). But, the fact 
that the counterpart is intrinsically bright, My ''^ —6.2 
would exclude the smallest stars. Looking again at Fig- 
ure H (right panel) wc sec that the SED of N5408 X-1 is 
well consistent with that of a BOI supergiant star (scaled 
at the magnitude of a lab class) and rules out all main- 
sequence stars. The turnoff in the F225W filter seems 
to match the spectral template of the BOI star. Also, a 
drop compared to a power law is expected in the NIR 
(Figure m right panel) and is clearly observed at l^m 
(Filter F105W). However, this drop should be even more 
prominent at longer wavelengths but is not really seen at 




Figure 4. Spectral energy distribution of NGC 5408 X-1 (corrected for extinction) for each observation (left), compared to stellar templates 
(right). SEDs 2 and 3 of NGC 5408 X-1 have been shifted down by 0.6 in log(Fi,) for clarity purposes. Fluxes of stellar templates in th e 
-ffST/WFC3 filters have been calculated using synphot and the Castelli and Kurucz stellar atmosphere models HCastelli fc Kurucdl200'3 ). 
In the right panel is ovcrplotted (black dot dashed line) an average power law fit from the three HST observations. 

1.5/iTO (Filter F160W). Proving that the companion star 
dominates the optical emission would be possible if opti- 
cal spectroscopy of the optical counterpart would reveal 
absorption lines, such as in the ULX P13 (jMotch et al.l 
120 lit) . But the signal to noise of th e optical spectra s tud- 
ied bv lKaaret fc Corbell (|2009t) and lCseh et all (I20T1I ) are 
not high enough to rule out the presence of an early B su- 
pergiant star because absorption lines in the wavelength 
range observed are not very deep. 



4.3. X-my to NIR SED 

The disk emission, extrapolated from two different 
models that we used to fit the X-ray data (a disk black- 
body with the addition of either a power-law or a thermal 
Comptonization model) does not fit the UV/optical/IR 
at all (Figure [5]) . For the power-law plus disk blackbody 
model, the extrapolated disk model underestimates the 
HST fiuxes by about an order of magnitude (Figure [5]) 
which basically means that the accretion disk emission 
would be negligible at optical wavelengths and the com- 
panion star would dominate. 

For the disk plus Comptonization model, the extrap- 
olated disk spectrum has a higher flux at lower energies 
but is still unable to explain the UV/optical/IR mea- 
surements (Figure [Sj . The excess emission could come 
from the donor star. However, the shape of the excess 
emission drops at the longest wavelengths, especially in 
the third observation where the NIR flux predicted by 
the accretion disk model is above the flux measured with 
HST (Figurc[5]). This is not consistent with the spectrum 
of any known star. 

Another possibility is that the emission in excess from 
the extrapolation of the disk may not come primarily 
from a stellar component. The irradiated disk model that 
we used is able to fit all the data quite well (see Table [3]). 
All three datasets give similar results within the errors, 
which is expected given the lack of significant variabil- 
ity both at X-rays and at lower energies. In this model, 
the UV/optical/IR emission is from the irradiated outer 
disk that thermalizes a fraction of the bolometric lumi- 
nosity foul ~ 0.03. This fraction is higher than seen from 
Galactic X-ray binaries in the thermal dominant state 



(jGierlihski et al.ll2009l: iZurita Heras et al.ll20lH) but sim- 
ilar fractions have been obtained from observations of 
the mi croquasar GRS1915 -I-105 in its most X-ray active 
p hase (iRahoui et al.l lMol). 

ICseh et al.l (I2011D have discussed the size of the ac- 
cretion disk in N5408X1, based on the broad emission 
lines seen in the optical spectrum of the counterpart. 
They concluded that the disk has a radius lower than 
2.35 X Mbh/1500 Mq AU which translates to 2.4 x 10"- 
10^^ cm for a 10-1000 M© black hole. According to the 
results from the irradiated disk model, we found that the 
disk would have a size '^l — 5x10^^ cm depending on 
the inclination of the disk. This is comparable to the size 
of the disk inferred in GRS1915-f 105 with - 3 x lO^^ ^m, 
where an irradiated accretion disk may also be present 
(jRahoui et al.ll2010[) . In N5408X1, the irradiated disk 
model gives i?out/^in ~ 1100-1600 which is i n the upper 
range of other ULXs, like in NGC 1313 X-1 (jYang et al.l 
[20T1I) with R^,,JR-, „ - 100-2000 and in NGC 6946 X-1 
( Kaaret et aL|[2010l ) with Ront/R\n - 40-6000. 

While disk irradiation is able to explain the observa- 
tional properties of N5408X1, we cannot make a defini- 
tive statement that irradiation is occurring in this sys- 
tem. First, no detection of significant, associated vari- 
ability between the UV/Optical/NIR and X-ray has 
been made (see below). Second, it is possible that the 
X-ray soft component is not emitted by an accretion 
disk, but by the photosphere of an outflow /wind (e.g., 
iPoutanen et al] 120071 : [Middleton et al.ll20lH) that would 
be due to supercritical accretion onto the accretion disk. 

In that context, it may be interesting to draw a com- 
parison with the X-ray transient V4641 Sgr. This is an 
interesting source because it may be a closer example to 
ULXs than low-mass X-ray binaries (LMXBs) thanks to 
its companion that is a quite massive star [M ^ 6.5 Mq). 
V4641 Sgr went into outburst in 1999 reaching an X- 
ray luminosity near or above the Eddington lim it of its 
- 10 Mq black hole (|Revnivtsev et all l2002al ) . Along 
with the X-ray outburst, the optical emission increased 
by at least 4.7 mag. It has been shown that irradiation 
in the accretion disk cannot be reconciliated with the ob- 
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Figure 5. Left: best fits for the simultaneous X-ray - UV/optical/NIR data of NGC 5408 X-1. The three plots correspond to the three 
independent observations. In each plot, the two accretion disk models (see text) extrapolated from the X-ray data fit to lower energies are 
shown in green (dot-dashed lino) and in blue (dashed lino), respectively corresponding to the disk model associated with a Comptonized 
component and the disk model associated with a power law component. They underestimate the UV/optical/NIR fluxes. The best fitting 
model is an irradiated disk (red line) that explains both the high energy and low energy data. This supports the interpretation that a 
standard, thin disk is present in NGC 5408 X-1. The spectrum of a BOI supergiant star (from the Castelli and Kurucz stellar atmosphere 
models, Castelli et al. 2004) is ovcrplotted in orange. It fits the HST data as well as the irradiated disk model. Note that the apparent gap 
between ~ 100 and ~ 1000 A is due to the extinction. Right: required flux needed to explain the HST measurements, after subtraction of 
the two extrapolated disk components seen in the left panels (the same color /linestyle is used to differentiate the two disk models). 

disk and thus would not contribute to the reprocessed 
emission in the optical bands. This is a natural con- 
sequence of the supercritical acc retion models because 
they predict geometric be aming (jPoutanen et al.l [20071 : 
lOhsuga fc Mineshigell20lll ). If the irradiated disk model 
is valid, it means that such beaming can be ruled out 
in N5408X1, because of the high fraction of thermalized 
bolometric luminosity in the disk. This would challenge 
the supercritical model. On the contrary, if the opti- 
cal spectrum is dominated by the donor emission, then 



served parameters. Instead, iRevnivtsev et al.l (l2002al lb[) 
argue that an extended envelope surrounding the source 
absorbs the X-ray flux and reemits it in the optical and 
UV. This would be the signature of a massive outflow 
driven by the supercritical accretion. The optical neb- 
ulae present around some ULXs and around N5408X1 
(Figure [p argue for the p resence of such outflows. As 
noted in lKaaret fc Corbell (|2009( ) the difhculty with the 
supercritical model is that the soft X-rays expected from 
the photosphere of the wind would not thermalize in the 
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Table 3 

Spectral Fit Parameters of the Irradiated Disk Model 



(kcV) 



(kcV) 



log (r 



"■^°-0.05 
, + 0.08 
-0.05 
fl + 0.07 
°-0.05 



"-470 



0.12;^ 
0.16 + 



2 42 + "- i« 

n + 0.17 
-0.19 

^o.l8 

^-0.19 



2.40_ 
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i, + II.U2 
"■^■'-0.02 
1+0.02 
-0.02 
-+0.02 
-0.02 



0.147 
0.14+ 



^■■^ -0.28 
, +0.66 
-0.37 
, + 0.25 
^-0.20 



0.713 
0.43 + 



"■"+^-0.022 

046 + '' °^'' 
O.U4b_Q P24 

" "'"'-0.018 



-0.25 
o ,,+0.28 
■'■^'--0.28 
7+0.27 
-0.29 



3.04^ 



3.07^ 



144.6 (93) 

101.7 (96) 
114.0 (103) 



Spectrum index used in the text 

Total absorption column, including the galactic extinction towards the source (njj = 0.057 X 10^^ cm~ 
Model normalization 
Power-law photon index 

Innermost temperature of the unilluminatcd disk 
Electron temperature 

Ratio of luminosity in the Compton tail to that of the unilluminatcd disk 
Radius of the Compton illuminated disk in terms of the inner disk radius 
Fraction of bolometric flux which is thermalized in the outer disk 
loglO of the outer disk radius in terms of the inner disk radius 



All errors are at the 90% confidence level. 

mild geometric beaming consistent with the constraints 
from the surrounding He II4686 nebula would be allowed. 
Thus confirming or not that the optical emission comes 
mainly from the donor star is important. 

We finally note that the reflection-based model 
used to model the X-ray sp ectrum of N5408X1 
(|Caballero-Garcia fc Fabianll2010[ ) would probably imply 
a negligible contribution of the accretion disk in the op- 
tical wavelengths. This is because there is apparently no 
need for thermal emission from the disk in that model. 
This may be explained by magnetic extraction of en- 
ergy from the disk, that would suggest that the energy 
is extracted before it is able to thermalize in the disk 
(|Caballero-Garcfa fc Fabianll201Cl[ )V In that context, ir- 
radiation in the outer parts of the disk would be at 
best limited, and the observed optical emission would 
be mainly due to the donor star. 

4.4. Variability 

Observing X-ray/optical correlated variability would 
be a direct test of the irradiated disk model. A lack 
of significant optical variability in response to a signif- 
icant X-ray variation would rule out the model, and 
would probably argue for the dominance of the com- 
panion star at non- X-ray wavelengths. The X-ray vari- 
ability we measure is limited, with a 40% change in the 
count rate between the first and last observation, and 
with the second observation displaying a flux consistent 
with that of the first one. At longer wavelengths, there 
is no clear, statistically significant variability in the HST 
filters. In practically all filters, the magnitudes in all 
three epochs are consistent within the errors (Table [2]), 
with only marginal evidence of variability between the 
first and third epoch, at a level < 2a. Thus, we can only 
provide an upper limit on the variability in the HST fil- 
ters, consistent with < 0.1 mag in most filters. If the 
situation is comparable as that of LMXBs, we would ex- 
pect the amplitude of optical variability to be roughly 
the square root of the amplitude of t he X-ray variabil- 
ity (jvan Paradiis fc McClintod^ 119941) . Our results are 
not inconsistent with this because the X-ray variability 
would imply a change of the optical flux of only ~18%, 
which is close enough to 10% and also because other pro- 
cesses can affect the optical luminosity of the system with 
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Figure 6. Absorbed X-ray flux (0.3-10 keV) versus normalization 
at 5500 A. X-ray fluxes are taken from the power law fits and 
normalizations come from the F336W/F547M/F845M fits. 

approximately the same amplitude (i.e., ellipsoidal vari- 
ations) . 

Perhaps the best evidence of variability may be seen in 
the normalization of the power-law fits. At 5500 A the 
difference in the normalization factors between the first 
and last observation is 0.67 ±0.23 x 10~^* erg s~^ cm~^ 
which is almost significant at a Scr level. We tried to 
look for a correlation between the X-ray flux and the 
magnitude/normalization in all fllters but the errors on 
both parameters are too large (see Figure [6|). 

The long term X-ray light curve of N54 08 X-1 mea - 
surcd with Swift shows dips quite regularly ()Konell201lD , 
with the X-ray luminosity decreasing by nearly an order 
of magnitude compared to the long-term average. In the 
context of an irradiated disk, this would imply a drop 
in the optical flux by a factor 3, i.e. ^ 1 m a g and 
thus would be easily measurable. iGierlinski et al.l (|2009f ) 
noted that this 'reprocessing signature' is approximate, 
depending on parameters such as the size of the disk, but 
given the order of magnitude of the variability that we 
expect, it should still be detectable, even if reduced, in 
the UV/Optical/NIR SED of the source. Finally, another 
additional signature would be an associated decrease in 
the flux of the He II emi ssion line visible in the optical 
spectrum of the source (ICseh et al.l[20ll . ii the emis- 
sion of He II is really associated with reprocessing in the 
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accretion disk. Thus, it would be of interest to perform 
optical observations during an X-ray dip. 

5. CONCLUSION 

We have studied NGC 5408 X-1 using 3 epochs of 
simultaneous Chandra/ HST observations. The optical 
counterpart of the ULX is visible from the UV to the 
NIR. The HST/WFC3 observations reveal that the souce 
is located near a young, ^ 5 Myr old, OB association 
containing massive stars up to 40 Mq, with a few blue 
supergiants. This could be the birthplace of the ULX 
system, but it cannot be ruled out that it was formed 
instead in a superstar cluster in NGC 5408. 

The UV to NIR SED is compatible with that of a 
BOI supergiant star. Using the simultaneous X-ray data, 
we show that the intrinsic emission of a standard accre- 
tion disk cannot explain the UV/Optical/NIR fluxes. In- 
stead, a model that takes into account irradiation in the 
disk fits well all the data. Further test of this model re- 
quires observing the source with a lower or higher X-ray 
flux and studying the correlation with the optical flux. 
This would be a good discriminant toward determining 
the physical origin of the optical emission. 
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Spectrum index used in the text 

External absorption column (in addition of the galactic extinction towards the source, nj] = 0.057 X 10^^ 
" Power-law photon index 
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All errors are at the 90% confidence level. 



